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ABSTRACT 


A  procedure  is  advocated  to  replace  the  common  practice  of  establishing  a 
rate  of  crack  growth  law  directly  from  increments  of  growth  or  from  slopes  of 
smooth  curves  drawn  through  crack  length  vs  number  of  loading  cycles  data.  In 
the  advocated  procedure,  a  postulated  rate  equation  such  as,  for  example,  the 
power  law  da/dN  'v  (AK)9,  is  integrated  for  the  specimen  geometry  used  and 
values  of  the  integral  corresponding  to  experimental  crack  lengths  are  plotted 
against  cycles  of  loading  to  reach  these  lengths.  Thus,  where,  and  if,  the 
rate  equation  is  valid  the  plot  will  consist  of  one  or  more  straight  line 
segments  for  each  of  which  constant  parameters  may  be  obtained,  or  confirmed 
if  theoretical.  Intersections  of  straight  line  segments  will  correspond  to 
crack  rate  of  growth  discontinuities  unobservable  by  the  usual  method.  Using 
Trip  Steel  data  supplied  by  Syracuse  University  as  an  example,  it  was  found 
that  all  growth  can  be  proportional  to  (AK)Z  with  a  proportionality  constant 
that  changes  discontinuously  at  various  amounts  of  growth  in  general  conform¬ 
ance  to  the  author's  experience  and  his  crack  growth  law.  This  conclusion  is 
not  invalidated  by  a  gradual  decrease  in  loading  during  testing,  or  by  the 
differences  in  analytical  expressions  for  AK  found  by  different  investigators. 
It  is  shown  that,  although  a  rate  proportional  to  (AK)1^  is  not  that  for  any  ' 
segment  of  the  whole  growth  curve,  it  represents  the  envelope  to  such  segments 
and  thus  is  a  sort  of  overall  representation  that  may  be  useful  in  design.  The 
discrepancy  between  a  rate  proportional  to  (AK)m,  m  =  3.7,  found  by  Syracuse 
University  and  the  rate  proportional  to  (AK)2  is  explained  by  showing  that  the 
Syracuse  method  of  analysis  gives  a  power  corresponding  to  an  envelope. 
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NOMENCLATURE 


CTS 

N 

B 

w 

a 

I  CL 
GCL 
x 

K 

P 

AP 

AK 

AP 

mean 

AK 

mean 

I 

gm 

ICTSm 

ICTSmD 

f(a/w) 


Abbreviation  meaning  compact  tension  specimen  of  ASTM  Standard  E399-70T 
Cycles  of  loading 
Thickness  of  CTS  (Figure  1) 

Width  of  CTS  (Figure  1) 

Crack  length  measure  for  CTS  in  ASTM  Standard  E399-70T  formula  for  K 
(distance  from  loading  axis  to  crack  front) 

Initial  (machined)  crack  length  of  CTS  (Figure  1) 

Grown  crack  length  of  CTS  (a-ICL)  (Figure  1) 

x  **  a/w  in  ASTM  CTS  formula  for  K  [Equation  (7)] 

Stress  intensity  factor 

Load  on  CTS 

P  -  P  .  in  CTS  loading  cycle 
max  min  J 

K  -  K  in  CTS  loading  cycle 
max  min 

(AP,  start  of  test  +  AP,  end  of  test)/2 
(AK,  start  of  test  +  AK,  end  of  test)/2 

Integral  for  growth  rate ~ (AK)m  for  geometry  g,  e.g.,  g  =  cc  for  a 
center  cracked  specimen 

CTS  integral  for  growth  rate~(AK)m 

I___  modified  for  decrease  (D)  in  load  during  test 
CTSm 

From  formula  for  AK,  i.e.,  AK  =  (AP/B-\/w)f  (a/w) 
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INTRODUCTION 


The  test  results  analyzed  here  were  supplied  by  D.  N.  Lai  and  V.  Weiss  of 
Syracuse  University,  who  obtained  them  to  satisfy  Navy- Syracuse  Contract  N00140- 
70C-0126. 

This  contract  was  concerned  with  the  mechanical  properties  of  Trip  Steels, 
including  toughness  and  crack  growth  characteristics. 

* 

In  this  connection  Lai  and  Weiss  made  an  experimental  study  of  crack  growth 
under  repeated  loading  using  an  ASTM  Compact  Tension  Specimen  (CTS).  They  analyzed 
the  growth  with  unusual  care  in  that  they  made  all  their  comparisons  of  the  rate  of 
crack  growth  for  a  specific  crack  length  of  0.6  in.  from  the  line-of-action  of  the 
loading  pins,  different  specimens  being  used  for  each  different  loading,  i.e., 
different  stress  intensity  parameter  AK.  On  each  test  the  minimum  load  used  was 
approximately  one-tenth  of  the  maximum  load. 

Their  main  conclusions  pertinent  to  the  present  study,  were  that  Trip  Steel 
offered  no  special  advantage  in  cyclic  crack  growth  resistance  and  that  the  rate 
of  crack  growth  was  almost  proportional  to  the  three  and  seven-tenths  (3.7)  power 
of  the  cyclic  range  of  toughness  parameter  AK. 

Not  only  was  AMMRC  interested  in  the  potential  of  Trip  Steels  in  general, 
but  in  formulae  for  the  toughness  parameter  K  and  the  crack  growth  characteristics 
in  particular.  Beeuwkes  of  AMMRC  had  personally  and  indirectly  analyzed  experi¬ 
mental  crack  growth  curves  of  a  large  number  of  materials  in  accordance  with  his 
crack  growth  formula  and  invariably  found  that  the  growth  could  be  very  closely 
matched  by  making  the  growth  rate  proporational  to  the  second  power  of  AK  in  con¬ 
stant  stress  amplitude  tests.  Thus,  it  was  of  interest  to  both  Syracuse  University 
and  AMMRC  to  determine  whether  or  not  Trip  Steels  behaved  differently  in  crack 
growth  caused  by  repeated  loading,  than  the  other  materials  analyzed  and,  in 
particular,  whether  the  rate  of  crack  growth  conformed  to  the  Beeuwkes  formula 
which,  essentially,  makes  yield  strength  the  only  adjustable  rate-controlling  mate¬ 
rial  property.  Accordingly,  Syracuse  University  turned  over  to  AMMRC  all  its  crack 
growth  data  on  this  steel  for  further  analysis.  This  was  carried  out  as  described 
in  this  report,  and  as  will  be  seen,  required  the  production  of  a  number  of  for¬ 
mulas  and  tables,  as  well  as  a  comparison  of  toughness  formulae  for  K  proposed  by 
different  sources,  all  of  which  may  be  useful  in  other  investigations  of  crack 
growth. 

It  is  found  that  the  results  obtained  often  strongly  depended  upon  the 
method  of  analysis  of  the  experimental  data.  Comparisons  are  made  with  the  pro¬ 
cedures  advocated  and  applied  here  to  all  the  test  results.  Accordingly,  this 
report  may  be  considered  to  be  a  primer  of  what  are  here  regarded  as  basic 


* 

The  Syracuse  findings  were  published  in  WEISS,  V.,  LAL,  D.N.,  and  BLOCK,  U. 
Mechanical  Behavior  of  Trip  Steel  as  a  Function  of  Fabrication  and  Processing 
Variables.  Dept,  of  Chem.  Engr.  &  Matls.  Sci.,  Syracuse  University,  New  York, 
MET-VW-1826-0272FR,  February  1972,  Submitted  to  Naval  Ship  and  Development  Lab., 
Annapolis,  Md.  21402,  Code  A933,  received  after  the  present  analysis  was  carried 
out. 
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concepts  of  a  proper  analysis  procedure.  Tables  to  facilitate  use  of  such  a  pro¬ 
cedure  are  provided. 


SYRACUSE  UNIVERSITY  TESTS  AND  METHOD  OF  ANALYSIS 

The  Syracuse  University  work  covered  the  mechanical  behavior  of  Trip  Steels 
as  a  function  of  fabrication  and  processing  variables.  Table  1. 


TABLE  1.  FABRICATION  AND  PROCESSING  OF  STEELS  TESTED  BY  LAL  AND  WEISS 


MATERIAL 

ALLOY 

CONDITION 

INITIAL 

TRIP 

DESIGNATION 

COMPOSITION 

PROCESSING 

PROCESSING 

Heat  2321 

C-0.27,  Mn-0.90, 

Final  Thickness 

Plate  hammer  forged 

80%  warm  roll 

Type  I 

Si-1.88,  Cr-8.80, 

0.300  in. 

1180C(2150F), 

426C(800F) , 

(B.) 

Ni-8. S,  Mo-4.0 

Austenitized  1204C 

Temper  350C 

1 

(2200F),  Cool  to  room 

(660F)  1  hr. 

temperature 

Yield  Strength:  220,000 

psi  longitudinal; 

200,000 

psi  transverse 

Heat  2322 

C-0.27,  Mn-0.91 

Final  Thickness 

Plate  hammer  forged 

80%  warm  roll 

Type  I 

Si-1.84,  Cr-8.81, 

0.28S  in. 

1180C(2150F) , 

426CC800F), 

(B  ) 

Ni-8. 73,  Mo-4 . 07 

Austenitized  1204C 

Temper  350C 

L 

(2200F) ,  Cool  to  room 

(660F)  1  hr. 

temperature 

Yield  Strength.  208,000 

psi  longitudinal; 

206,000 

psi  transverse 

Heat  2322 

C-0.27,  Mn-0.91, 

Final  Thickness 

Plate  hammer  forged 

80%  warm  roll 

Type  II 

Si-1.84,  Cr-8.81, 

0.31S  in. 

1180C(2150F) , 

426C(800F) 

(B.) 

Ni-8. 73,  Mo-4.07 

Austenitized  1204C 

15%  Cold  work, 

J 

(2200F) ,  Cool  to  room 

Temper 

temperature 

3S0C(660F) ,  2  hrs. 

Yield  Strength:  220,000  psi  longitudinal; 

208,000  psi  transverse 

.  . . 6 

Note:  In  the  present  report.  Young's  Modulus  of  Elasticity  was  assumed  to  be  E*28  x  106  psi. 


The  nominal  composition  was  0.27  C,  0.90  Mn,  1.86  Si,  8.80  Cr,  8.7  Ni,  and 
4.0  Mo.  Two  heats  were  employed  and  two  Trip  Steel  processes  were  used  on  one  of 
the  heats  giving  a  total  of  three  materials: 


1. 

Bl 

(Heat 

2321, 

Type 

I) 

2. 

b2 

(Heat 

2322, 

Type 

I) 

3. 

b3 

(Heat 

2322, 

Type 

II). 

All  three  materials  were  tested  for  crack  growth  characteristics.  As  a  result 
it  was  concluded  in  the  Syracuse  University  work  that  a  fatigue  crack  propagated 
through  the  material  according  to  the  rather  generally-used  power  law  rate  of  crack 
growth  per  cycle  relationship. 


da 

dN 


(1) 


-  C  (AK)m 
m 


with  m  =  3.7.  Here  AK  is  the  maximum  value  of  K  minus  the  minimum  value  of  K  in 
the  loading  cycle.  (The  Beeuwkes  relationship  that  we  shall  employ  in  part  has 
m  *  2  and  a  theoretical  value  for  C.) 

Lai  and  Weiss'  test  method  consisted  of  fatiguing  about  four  ASTM-types  of 
CTS  (Figure  1)  of  each  material  at  different  stress  intensity  amplitudes  AK  with 
R  3  minimum  load/maximum  load  =  0.1  until  the  fatigue  crack  had  grown  0.4  in. 

Figure  2  shows  four  samples  of  crack  growth  curves  and  data  supplied  to  the  writer. 
Table  2  presents  data  on  the  circled  experimental  points  from  which  the  crack 
growth  curves  were  drawn  and  Table  8  presents  the  data  from  these  curves. 

Since  the  load  decreased  somewhat  during  testing,  load  range  values  were 
computed  for  the  beginning  and  end  of  the  entire  crack  growth  curve.  These  two 
values  were  averaged  for  computation  of  the  stress  intensity  amplitude  AK  used  in 
analysis  of  the  data.  The  da/dN  values  used  were  determined  from  tangents  to  the 
curves  at  a  fixed  "grown"  crack  length  of  0.2  in.,  i.e.,  "a"  in  Figure  1  was 
0.6  in. 

Their  resultant  plot  of  AK  vs  da/dN  is  shown  in  Figure  3  in  comparison  with 
other  work  they  referenced. 

PRESENT  TYPE  OF  ANALYSIS 

General  Considerations  and  Power  Crack  Growth  Law 

A  postulated  rate  of  crack  growth  equation  is  integrated  here  to  give  a 
straight  line  relationship  between  the  integral  and  the  number  of  cycles  occuring 
between  the  limits  of  integration.  If  the  form  of  rate  equation  is  valid,  values 
of  the  integral  corresponding  to  observed  total  crack  lengths  plotted  against  the 
loading  cycles  to  reach  these  crack  lengths  will  lie  on  a  straight  line  whose 
slope  may  be  compared  with  a  theoretical  one.  It  is  felt  that  this  method  is 
superior  to  the  usual  approach,  that  of  trying  to  correlate  with  the  unintegrated 
da/dN  rate  equation  itself,  largely  because  of  difficulties  in  ascertaining  satis¬ 
factory  experimental  values  of  da/dN.  If  da/dN  is  taken  to  be  Aa/AN  where  Aa  is  a 
small  increment  of  crack  growth,  it  is  difficult  to  observe  a  suitable  Aa  for  small 
AN  because  of  irregularity  and  lack  of  definition  of  the  crack  front.  In  fact,  Aa 
may  be  different  according  to  its  observed  location  (sides  of  test  piece  or  inter¬ 
ior).  Furthermore,  it  is  difficult  to  obtain  an  appropriate  rate  from  the  slope  of 
a  curve  drawn  through  selected  da/dN  data  because  of  the  diverse  characteristics 
and,  especially,  lack  of  possible  abrupt  changes  in  slope,  that  apparently  equally- 
suitable  matching  curves  may  have.  In  the  present  method  only  straight  lines  are 
run  through  the  data;  either  the  integrated  rate  equation  is  suitable  for  this  or 
it  must  be  changed.  As  will  be  shown,  the  rate  at  any  crack  length  is  very  simply 
given  by  the  one  slope  of  the  line  through  the  range  of  matching  data  together  with 
the  value  of  AK  at  that  crack  length. 
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DIMENSIONS  GIVEN  IN  INCHES 


0.30(B,) 
0.285(B2) 
0.31 5(B3) 


Lal-Weiss  reportedly  used  the  A  STM  standard  E  399-70T  CTS-  geometry  for  plane  strain  fracture 
testing  of  metallic  materials.  However,  some  of  their  dimensions  substantially  differ  from  this,  as 
shown  below,  omitting  tolerances.  ICL  varied  somewhat  and  is  noted  on  the  data  for  each  specimen 
All  ASTM  standard  distances  are  proportional  to  w.  "a"  is  the  ASTM  standard  crack  length.  ICL 
is  the  initial  crack  length. 


Dimension 


Pin  Hole  Diameter 
Height  ( A  Crack) 
Width 
Thickness 
Pin-to-Pin 
Slot  to  Top 


Lal-Weiss 


0.375  in. 
2  in. 

2.20  in. 
~0.3  in. 
1  in. 

1  in. 


ASTM  standard  E  399-70T 


0.25W  =  0.45  in. 
1.2W  =  2.16  in. 
1.25W  =  2.25  in. 
0.5W  =  0.9  in. 
0.55W  =  0.99  in. 
0.6W  =  1.08  in. 


MATERIAL  -  B,(^{AT  2321, TYPE  I  ) 
SPECIMEN  No.  I 


MATERIAL  -  B2  (HEAT  2322  .TYPE  I  ) 
SPECIMEN  No. 3 


AT  GROWN  CRACK  LENGTH  OF  0.2  In. 


A  K  =  21.2  KSI  j/in. 

=  5.8  h-  in./CYCLE 


ATGROWN  CRACK  LENGTH  OF  0.2  in. 


O  20,000  40,000  60,000  80,000  100,000 
NUMBER  OF  CYCLES 
Pmax  (initial)  =  2,000  lb  1  A 
Pmox  (final )  =  1 ,000  lb  j  A  ” =  °p  P™'{,n*on)  * 


Pmox  (final )  =  1 ,000  ,b)  A  P  =A°p9  = '.™  »b 

pmin  =  0.1  Pmax  AK*b7t,(*) 

B=0.300in.\  , 

W=l.800in.)fl")~5 
INITIAL  CRACK  LENGTH  =  0.394  in. 

(a) 

MATERIAL  -  B3  (HEAT  2322, TYPE  H  ) 

SPECIMEN  No. 4 

g'm~-  AT  GROWN  CRACK  LENGTH  OF  0.2  in. 


-0.40 

x 


AK  =  30.6 KSI  t/in. 

=  20/.  in  /CYCLE 


£  0.40 

X 


w  0.30 

x 

o 

g  0.20 

z 

» 

o 

x  0,10 
o 


AK=  27.2  KSI i/inT 

■■f  =  16.3/1  in./CYCLE 
ON 


0  5,000  10,000  15,000  20,000  30,000  35,000 

NUMBER  OF  CYCLES 

Pmox  (initial): 2,500 lb  AP(mean)  =2,080  lb 
Pmox  (final  )=  2,122  lb  8=  0.285  in 
Prmn=O.IPmax  W=  1.000  in. 

INITIAL  CRACK  LENGTH*  0  400  in. 


MATERIAL  -  B3  (HEAT  2322,  TYPE  H  ) 
SPECIMEN  No.  5 


AT  GROWN  CRACK  LENGTH  OF  0  2in. 


AK=  32.6  KSI  i/in7 
^  =29/iin. /CYCLE 


•S  0.40 


X 

o 

< 

o  0.20 


9 

o 

S  o.io 


uj  0.30 


x 

o 

2  0.20 
o 


I 

g  O  iO 


0  5,000  10,000  15,000  20,000  25,000 

NUMBER  OF  CYCLES 

pmax  (initial )  =  3,0001b  AP(mean)  =  2, 589.5  lb 
PmaX(final)  =  2,7541b  B=0.3l5in. 

pmin  :  0.1  Pmox  W  =  l.800ln. 

INITIAL  SLOT  LENGTH  =  0.400 in. 


0  5,000  10,000  15,000  20,000  25,000  30,000 

NUMBER  OF  CYCLES 

Pmox(initiol)  -  3,200  lb  AP(mtar)  ;  2,7581b. 
Pmax(finol)  =  2,9291b  B=03l5in. 

Pmin:O.IPmox  W=l.800in. 

INITIAL  SLOT  L  E  NGTH  =  0.405  in. 


Figure  2.  Four  samples  of  data  sheets  supplied  by  Lal-Weiss  (3/8  X  actual  size). 


ORIGINAL  TRIP  STEEL  DATA  OF  GROWN  CRACK  LENGTH  (GCL)  vs  CYCLES  (N)  AND  THE  CORRESPONDING  m  =  2  CTS  INTEGRAL  I 

USING  ASTM  STANDARD  E399-70T  FORMULA  FOR  K 


Figure  3. 


I 

Log-log  plot  of  fatigue  crack  growth  rates  da/dN,  vs  stress  intensity  amplitudes, 
AK,  of  Lal-Weiss  and  Gerberich.1 


Thus,  consider  Equation  (1) 


da 

dN 


C  (AK) 
m 


m 


(1) 


in  which  the  constant 


t 


Cjj,  evidently  has  the  dimensions  (levgth) 


l-m/2 


(stress) 


Equation  (1)  may  also  be  written  as 


d (a/b)  I 

(  AK  \ 

m 

C  bm/2(AS  )m 
m  n 

dN  1 

.l/2._ 

ASn/ 

b 

(la) 


Note  that  the  stress  dimension  in  Cjj,  is  not  determined  when  Cm  is  from  tests  made 
with  a  single  value  of  ASn.  For  example,  unless  otherwise  determined,  the  stress 
in  Cfl,  might  be  assumed  to  be  [exYY (ASn) l-x-y]  where  E  is  the  Modulus  of  elasticity 
and  Y  is  yield  strength;  thus,  in  this  case  da/dN  is  not  really  proportional  to  AK 


1.  GERBERICH,  W.W.  Metastable  Austenitic  Steel  With  Ultra  High  Strength  and  Tough¬ 
ness.  SAE  Paper  No.  690262,  Int.  Auto.  Eng.  Cong.,  Detroit,  Mich.,  January,  1969. 
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where  b  is  an  arbitrary  length  dimension  that  factors  out  along  with  the  nominal 
loading  stress  range  ASn,  from  the  analytical  expression  for  AK  for  the  test  spec¬ 
imen  geometry  and  loading.  Thus,  AK/(b''2  ASn)  is  a  dimensionless  function  of 
(a/b)  and  dimensionless  ratio  constants  that  represent  the  test  specimen  geometry. 
(See,  in  particular,  ASTM  Standard  E399-70T  where  we  chose  b  *  w.) 

Hence, 


(a/b)  N 

1  I  f  J - d(°/b)  T .  -  b^2-1  C  f  US  )'  1>  E  l-'2'1  C  £ 

ci/b)  K<bl/2 *v]  1 


so  that  if  AS  is  constant 
n 


I  =  bm/2-1  C  (AS  )m(N  -  N  ). 
gm  m  n  o 


Thus,  if  C  is  constant,  Igm  is  rectilinear ly  ("straight  line")  related  to 
Em  or,  more  simply,  if  ASn  is  also  constant,  to  N. 

We  characterize  here  the  integral  on  the  right  by  Em  and  the  integral  on  the 
left  by  Igm,  where  g  stands  for  the  geometry  of  the  test  specimen  and  m  for  the 
power  of  the  rate  relation.  Equation  (1).  Thus,  Ic_s2  is  the  integral  for  the  CTS 
with  m  =  2.  The  lower  limit  (a/b)  can  have  any  value  so  long  as  it  is  below  the 
experimental  range  of  (a/b).  We  hive  provided  tables  of  IrTo->  and  Iptca  with 
(a/b)Q  =  xD  -  0.2  for  use  with  the  CTS. 

For  a  very  simple  example  of  I  consider  the  case  of  a  central  crack  (g  =  cc) 
in  a  wide  plate  under  uniform  tensile  load  perpendicular  to  the  crack  and  far  from 
it,  with  m  =  2.  Then, 


AK  =  ASVira 


so  that 


■  {  7^77 TV  '  <ln  a/a»,/u 


Here  we  have  chosen  b 


a  to  be  the  initial  crack  length, 
o 
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Note  that 


[1/(N 


l‘l/m 


(AS  )“  dN 
n 


(4) 


is  a  constant  stress  range  equivalent  to  the  actual  spectrum  AS  in  N  -  N  cycles 
of  loading.  n  ° 

It  was  evident  from  Equation  (2)  that  I  is  rectilinear ly  related  to  I,  or 
more  simply  to  N  where  ASn  is  constant.  Suppose  values  of  Igm  corresponding  to 
experimental  values  of  x  *  a/b  were  plotted  against  values  of  Z,  or  N,  correspond¬ 
ing  to  a/b  and  found  to  conform  to  a  straight  line  over  some  range.  Then,  in  that 
range. 


C  = 
m 


(I  )  -  (I  ) 

gm  2  gm  1 

(I  )  -  (z  ) 

m2  m  1 


,m/2-l 


(5) 


or,  where  AS  is  constant , 
n 


^Tgm^2  (Ig®\ 


m 


N2  "  N1  bm/2_1(AS  )m 


=  Slope^b”/2  1  (ASn)m]  , 


(5a) 


where  No,  (I  )  and  Ni,  (I  )  are  any  two  points  on  the  line  and 

gm  2  g®  1 

f(I  )  -  (I  )  1/(N2  -  Ni)  is  its  slope.  Thus,  C  is  determined  from  the  data  and 

l  gm  2  g®  lj  m 

we  see  that  a  straight  line  plot  of  I  vs  N  along  with  Equation  (5)  or  (5a)  com¬ 
pletes  the  match  to  Equation  (1). 

From  this  line  the  rate  of  growth  for  any  value  of  AK  is  then  found  by  sub¬ 
stituting  C  from  Equation  (5)  or  (5a)  into  Equation  (1),  e.g., 


-  b (Slope)  (ak/Oj1/2  AS  )lm 


(6) 


where  ASn  is  constant  or  the  constant  equivalent  to  the  actual  spectrum  in  the 
region  of  constant  slope. 

Experience  has  shown  that  I  vs  N  curves  are  often  made  up  of  straight  line 
segments  with  different  slopes.  Thus,  In  da/dN  vs  In  AK  curves  will,  in  accord¬ 
ance  with  Equation  (6),  be  made  up  of  straight  line  segments  with  a  jump  in  da/dN 
at  the  end  of  each  segment.  There  appears  to  be  a  tendency  for  these  transitions 
to  be  spaced  at  equal  intervals  of  In  AK  so  that  the  ratio  of  rates  at  the  transi 
tions  as  well  as  of  the  AK's  at  two  successive  transitions  is  constant.  In  this 
case  the  transition  points  are  connected  by  two  parallel  straight  line  envelopes 
(Figure  4).  There  also  seems  to  be  a  tendency  for  the  segments  to  lie  between 
these  two  envelopes  even  where  the  transitions  do  not  occur  at  exactly  equal 
intervals  of  AK. 


Figure  4.  Idealized  da/dN  vs  K  Diagram. 


In  Figure  4  let  da/dN  on  the  right-hand  envelope  be  r^  and  along  the  left- 
hand  envelope  be  r Then, 


where  p  and  the  8's  are  constant. 

Let  subscript  i  denote  a  specific  transition  and  let  C  .  be  C  of  Equation 
(1),  between  i  and  i  -  1.  Then,  *  H1 


rR, i  *  Cl-l,i<AK)"  ‘  *R<“>i 
rL,i-l  -  Cl-l,i(4K)"-l  '  *L(4K)Pi-l 


(7a) 

(7b) 


so  that  the  envelope  on  the  right  is 


rR  =  Vi.i^rw 


(8a) 


and  the  envelope  on  the  left  is 


rL  "  Ci-l,i(AK)^Pl(AK)P  • 


(8b) 


From  Equations  (7a)  and  (7b)  the  ratio  of  B's  as  well  as  the  ratio  of  rates  at 
the  end  and  beginning  of  a  segment  are  constant  since  the  ratios  of  the  correspond¬ 
ing  (AK)'s  are  constant,  i.e.. 


(AKV 


(AK) 


i-1 


m-p 


and 


R,i 


L,  i-1 


(AK) , 


(AK) 


i-1 


m 


(9) 


At  1,  (AK)  *  (AK) ^  so  that  from  Equations  (8a)  and  (8b)  the  jump  rate  ratio  is 


L,i+1 


R,i 


(AK). 


(AK) 


i-1 


p-m 


(10) 


This  is  also  the  ratio  of  the  C's  on  two  successive  segments,  i.e., 


Finally,  the  ratio  of  the  jump  rate  ratio  to  the  ratio  of  rates  at  the  end 
and  beginning  of  the  segment  ending  at  the  jump  is 


(rl.i+l/rR,l)/<ri/ri-l> 


(AK)i 


p-2m 


(12) 


which  is  a  constant. 

Analyses  of  many  experiments  suggest  that  m  »  2  and  p  =  4  and  that  the  ratio 
of  rates  at  the  end  and  beginning  of  any  m  =  2  straight  line  segment  spanning  the 
distance  between  the  envelope  lines  of  a  plot  of  log  da/dN  vs  log  (AK)  is  two  as 
specified  in  the  generalized  Beeuwkes  crack  growth  law  discussed  in  the  next 
section.  For  such  a  case  we  consider  a  segment  whose  equation  is  known  from  theory 
or  experiment  and  which  begins  at  a  known  value  (AK)  =  (AK)^^,  i.e.. 


d»  ■  ci-i,t("°2  ”lth 


(AK)^!  <  (AK)  <  (AK)± 


(13) 


where  ^  and  (AK)^_  are  known.  In  this  case  we  will  now  find  the  equations  of 

the  envelopes  and  the  characteristics  of  the  actual  growth. 

From  Equations  (7a)  and  (8a) 


rR,i 

rL,i-l 


(AK)  t 


(14) 


so  that 


(AK)±  -  J2  (AK)i_1 


(15) 


Hence  from  Equations  (8a)  and  (8b)  the  rate  equations  of  the  envelopes  are 


rL  ’  Cl-l,i  WK)4 


rR  ■  <“>1-1  <4K)4/2  ' 


(16a) 

(16b) 


so  that  log  r  -  log  r 


log  2. 


(17) 


The  ratio  of  the  rate  constants  of  two  adjacent  growth  segments  Is 


Ci+l,i+2  rL>i+l  _ 
Ci-l,i  rR,i  SR 


(18) 


The  ratio  of  the  jump  in  rate  from  one  segment  to  the  next  to  the  increase  in 
rate  from  the  beginning  to  the  end  of  the  first  segment  is 


<rL, i+l/rR, iJ / *ri^ri-l* 


4-2(2) 


Beeuwkes  Crack  Growth  Law 

This  crack  growth  law  for  applied  stress  amplitudes  (well  below  the  yield 
strength)  and  for  plane  strain,  the  present  case,  was  originally  written 


da 

dN 


(§)’(?) 


(20) 


and  may  be  written 


(20a) 


where 


and  where  u  is  Poisson's  ratio,  E  is  the  Young's  Modulus  of  Elasticity,  Y  is  the 
yield  strength,  and  S  is  the  "nominal  field"  stress,  i.e.,  the  stress  equivalent 
to  the  nominal  loading  stress  in  a  long,  wide  tension  specimen  having  a  small, 
central  crack  of  length  2a  lying  perpendicular  to  the  loading  direction,  i.e., 

S  *  Kj/tt  in  modern  stress  intensity  factor  notation.  In  the  case  of  considerable 


For  sharp  notches  other  than  the  small  central  crack  across  the  tension  field,  S, 
inawide  plate,  S  «  [kg/(2  ,/a/p)]  Sgn  as  p  +  0,  where  kg  is  the  elastic  stress 
concentration  factor  for  the  geometry  g  and  Sgn  is  the  corresponding  nominal  stress 
for  that  geometry  (k  includes  an  a/p  factor). 

*f*  ^  _ _ 

Taking  as  definition,  K  ■  S  ,/ri  in  the  case  of  the  above  tension  member  with  the 
central  crack. 
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plasticity,  S/E  was  to  be  regarded  as  a  loading  strain.  The  effect  of  different 
ranges  of  stress  had  not  been  worked  out. 


Subsequent  analyses  of  experiments  found  in  the  literature  indicated  that: 

(1)  at  least  for  tension-tension  stressing  (S^  -  Smin)  5  ASn  could  be  substituted 

for  S  so  that  (AK)2/tt  could  be  substituted  for  S^a;  (2)  that  at  the  start  of  most 
of  the  tests  the  constant  da/dN/(AK)2  was  uniformly  distributed  within  1/  -Jl  to 
times  Cg  in  Equation  (20a);  (3)  that  though  the  data  could  be  well-matched  assum¬ 
ing  the  rate  remained  proportional  to  (AK)2  for  constant  (Smax  -  Smin),  the  con¬ 
stant  of  proportionality,  Cg,  changed  abruptly  after  various  amounts  of  crack 
growth,  sometimes  by  a  slight  amount,  but  where  the  change  was  substantial,*  by 
factors  averaging  close  to  two  and  this  change  was  not  typically  associated  with  a 
transition  from  flat  to  shear  mode;  and,  finally,  (4)  that  though  the  length  of 
crack  at  changes  was  variable  even  under  similar  conditions  of  testing,  there  was  a 
tendency,  in  plots  of  log  da/dN  vs  log  AK,  for  the  second-power  straight-line  seg¬ 
ments  matched  with  the  substantially  different  constants  of  proportionality,  to  lie 
within  two- fourth-power  straight  line  envelopes. 

In  the  latter  case  the  ratio 


(AK  )2/( AK  )2  =  2 
ne  no 


where  ne  corresponds  to  the  end  and  no  to  the  beginning  of  segment  "n".  With  AS 
constant  this  is  (a^  /anQ)  =  2,  the  ratio  of  crack  lengths  at  the  beginning  and  end 
of  any  segment,  andnfhengrowth  is  thus  in  conformance  with  the  formulae  at  the  end 
of  the  preceding  section.  Equations  (13)  to  (19). 

For  short  cracks  the  uncertainties  in  crack  length  determination  are  sufficient 
to  make  the  length  of  the  line  segments  (a  -  a  )  =  a  indistinguishable  so  that 
the  experimental  law  of  growth  appears  to  61  tha6°of  the  fourth  power  of  AK  coupled 
with  (i.e.,  data  within)  a  scatter  band.  For  longer  crack  lengths  however,  cases 
have  come  to  the  present  writer's  attention  where  a  law  of  growth  has  been  suggested 
on  the  basis  of  experimental  data  that  seems  to  him  to  actually  correspond  to 
enclosing  envelopes,  the  experimental  points  of  individual  tests  making  up  the 
whole  plot  appearing  to  correspond  to  a  lesser  slope  than  that  of  such  envelopes. 

Recapitualting,  with  AS  =  S  -  S  ,  where  S  represents  nominal  applied 
loading  stress,  we  would  expect 


Higher  rates  are  natural  if  cleavage  areas  occur  at  the  AK  used,  as  has  been  ob¬ 
served  in  some  aluminum  tests,  but  there  are  other  possible  mechanisms  for  the  rate 
transitions.  For  example,  the  writer  has  suggested  that  growth  may  be  concentrated 
in  statistically  uniformly  distributed  points  along  the  crack  front  so  that  at  a 
transition  to  a  higher  rate  a  growth  point  in  between  each  is  added,  thus  doubling 
the  rate. 
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where  the  constant  Cg  introduced  here  is  expected  to  have  values  given  by 


2n/  y/l  <  CB/  [(2/3)1/3(l  -  p2)4/3]  <  2ny2  (22) 

instead  of  simply* 

CB  -  (2/3)1/3(l  -  y2)4/3  Sir /4  , 

the  value  of  C  in  Equation  (20a) . 

B 

Here  n  ■  0,  basically  but  there  may  be  transitions  with  n  having  the  values 
•  • • i  —  2 ,  *1|  0j  1,  2|  • » •  •  (23) 

If  n  «  0,  we  have  approximately 

it/ (4  Jl )  S  0.555  <  CB  <  1.111  =  ir^A  .  (24) 

In  Equation  (21),  the  term  (AK/AS  )  is  the  analytical  expression  for  AK  with  the 
nominal  stress  factored  out.  11 

SECOND  POWER  CTS  INTEGRAL  AND  ANALYSIS;  AK  BY  ASTM  STANDARD  E399-70T 
Constant  Loading  Range 

AK  for  a  CTS  can  be  calculated  from  the  power  series  law  given  in  ASTM  Standard 
E399-70T.  It  is 

AK  -  (AP/BWl/2)(x1/2)(29.6  -  185. 5x  +  655. 7x2  -  1017x3  +  638. 9x4)  ,  (25) 


n  was  chosen  to  put  scatter  in  experimental  determination  of  this  basic  value  C 
to  either  side  of  CB.  If  n  were  replaced  by  (n  +  1/2)  the  scatter  would  be  to  one 
side: 

2n  i  CB/  [(2/3)  1/3 (1  -  u2)4/3]  <  2n+1 
so  that  for  the  basic  form,  n  -  0, 

1  <  CB/  [(2/3)l/3(l  -  112)4/3]  <  2  , 
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where  with  a  -  distance  from  the  loading  axis  to  the  crack  front  (i.e., 
a  -  ICL  +  grown  crack  length).  Figure  1, 


x  *  a/w 


AP/BW  «  AS 


The  second  power  CTS  integral  expression  corresponding  to  this  formula  for  AK 
may  be  gotten  from  Equation  (2),  i.e.. 


x(29. 6  -  185. 5x  +  655. 7x2  -  1017x3  +  638. 9x4)2 


The  integral  was  calculated  and  tabulated  along  with  Iqts2  ^ rom  other  proposed 
expressions  for  K  to  be  discussed  later,  page  30,  Table  7,  briefly  in  Table  3  and 
more  fully  in  Appendix  A  for  the  range 

0.2  <  x  5  a/w  £  0.8 

The  lower  bounding  limit  on  x  was  selected  because,  for  CTS,  K  for  cracks  for 
x  =  a/w  <  0.2  is  too  sensitive  to  the  actual  distribution  of  applied  load  to  warrant 
general  formulation.  (Thus,  experimental  determination  of  K  for  x  <  0.2  for  each 
particular  loading  pin  configuration  being  used  is  desirable.)  The  upper  limit  was 
chosen  to  avoid  getting  into  a  region  of  side  edge  effects.  The  integral  is  plotted 
in  Figure  5. 


It  is 


Also,  a  Pade  approximant2  was  approximately  fitted  to  the  tabulated  results. 


1  -  3.400x  ~  8.030X2 
1  -  1. 364x  +  7. 04 Ox2 


X  5.548  X  10 


Using  Appendix  A  to  get  from  a/w  values,  vs  N  plots  were  made 

using  the  actual  "grown"  crack  length  vs  cycle  N  data  points  (Table  2),  as  well  as 
data  points  taken  at  regular  intervals  from  smooth  curves  drawn  by  Lal-Weiss 
through  the  actual  data  points  (Figure  2). 


BAKER,  G.  A.  JR.  Accurate  Long-Range  Extrapolation  -  The  Pade  Approximant. 


Brookhaven  National  Laboratory,  Brookhaven  Lecture  Series,  BNL  50241  (T576), 
Associated  Universities,  Inc.,  Contract  with  U.S.  Atomic  Energy  Commission. 


16 


Table  3.  CTS  INTEGRAL  I 
K  vs  X  s  a/w.  SECOND- 


tGKAL  iCTS2  x  TO  vs  x  =  a/w  F0R  AS™  STANDARD  E399-70T,  SRAWLEY,  BOWIE,  and  BEEUWKES 
SECOND-DEGREE  PADE  USED  FOR  BEEUWKES  AND  BOWIE  K,  THIRD  DEGREE  PADE  FOR  SRAWLEY  K. 
[SEE  EQUATION  (27),  FIGURE  5,  AND  TABLE  7.]  CONDENSED  TABLE. 


STANDARDS  SRAWLEY 
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It  was  found  that  the  IcxS2*  ^  points  were  lined  up  on  one  or  more  straight 
line  segments3  having  different  slopes,  even  for  data  taken  from  the  smooth  curves 
where  the  transitions  in  slope  corresponding  to  the  change  from  one  straight  line 
segment  to  another  was  not  apparent. 

Despite  the  scatter  in  the  actual  data  there  was  no  substantial  difference  in 
the  straight  line  representations  made  with  the  two  data  sources.  Thus,  in  the 
interest  of  brevity  and  economy  only  those  plots  made  with  the  actual  data  points 
are  reproduced  in  Figures  6  through  18. 

Since  the  plots  are  made  up  of  one  or  more  straight  line  segments,  the  crack 
growth  is  correspondingly  represented  as  increasing  as  the  square  of  the  crack  tip 
stress  intensity  range,  AK,  with  a  proportionality  constant  that  changes  abruptly, 
but  only  when  growth  shifts  from  any  segment  to  another. 

The  constants  (previous  section)  derived  from  Figures  6  through  18,  as  well 
as  those  derived  from  the  smooth  curve  data  not  reproduced  here  are  shown  in  the 
accompanying  Table  4.  Considering  irregularities  in  the  data  and  other  factors  to 
be  discussed  later,  the  agreement  with  present  type  of  analysis  is  considered 
satisfactory. 


3.  TREMBLAY,  R.J.  Fatigue  Crack  Growth  Transitional  Behavior.  Army  Materials  and 


Mechanics  Research  Center,  AMMRC  TN  73-13,  October  1973.  For  copies  of  this 
Technical  Note,  contact  Dr.  Beeuwkes,  AMMRC,  923-5744. 
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a/w 

Figure  5.  CTS  Integral  l^-ye  vs  a/w  based  on  ASTM  Standard 
E399-70T  formula  for  K,  and  constant  loading  range. 


Effect  of  Gradual  Decrease  in  Load  During  Testing 

The  matching  formulae  used  to  this  point  are  designed  for  constant  maximum  and 
minimum  loads.  Actually,  these  loads  decreased  to  around  seven-eighths  of  their 
initial  values  during  the  testing,  and  in  applying  the  previous  analyses,  the  mean 
of  the  initial  and  final  values  were  used.  (This  mean  was  also  used  by  Syracuse 
University  in  its  analysis.)  The  error  introduced  by  this  approximation  was  evalu¬ 
ated  as  follows,  and  so  found  to  be  essentially  negligible.  The  actual  load  range* 
was  taken  to  be 

AP  -  AP  (-0.6  a/w  +  1.198)  (29) 

mean 


with 


0.22  <  a/w  <  0.44 


Despite  considerable  scatter.  (See  r  in  Table  4.) 


25,000  50,000  75,000  100,000 
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Figure  6.  IcjS2  vs  N  for  material  B-j,  Type  I,  Specimen  1. 


B,  (2321,1) -2  2/2/73  LC 
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Figure  7.  IqtS2  vs  ^  for  mater'al  B,,  Type  I.  Specimen  2. 
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Figure  8.  1 0752  vs  ^  *or  mater>al  B-|<  Type  I,  Specimen  3. 
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Figure  9.  vs  n  f°r  material  B1#  Type  I,  Specimen  4. 
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ial  B2,  Type  I,  Specimen  1. 


B2  (2322,I)-2 
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B2,  Type  I,  Specimen  2. 
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N  for  material  82,  Type  I,  Specimen  3. 


B2  (2322, 1)-4 


2/2/73  LC 


N  CYCLES 


Figure  14. 


ICTS2  vs  ^  *or  mater'a*  B3.  Type  II,  Specimen  1. 
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vs  N  for  material  Bg,  Type  II,  Specimen  3. 
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Material 

and 

Specimen 

Max. 

Load 

Mean 

Range 

AP.lb. 

Fall-Off 

Ratio 

r 

Slope 

No. 

Original  Data 

Points 

Data  Sheet  Curves 

CB 

Slope  XI 0^  Growth 

Figures  6-18  Constant  Cg 

Initial 
Pxi *  1 b - 

Final 

P  f,lb. 

xf 

B !  I  -  1 

2000 

1800 

1710 

0.10 

1 

7.46 

0 

.7242 

0.4737,0.8969,0.5892 

BVI  -  2 

1500 

1300 

1260 

0.13 

1 

3.56 

0 

.6367 

0.6544 

2 

1.94 

0.3469 

0.4094 

B !  I  -  3 

2500 

2375 

2192 

0.05 

1 

6.66 

0 

.3934 

0.7207 

2 

17.07 

1 

.0082 

1.1165 

3 

24.72 

1 

.4602 

1.4946 

Bjl  -  4 

3000 

2782 

2602 

0.07 

1 

9.96 

0 

.4176 

0.4318 

2 

24.92 

1 

.0448 

1.0733 

B2I  -  1 

2000 

1566 

1685 

0.22 

1 

18.38 

1 

.5977 

0.5728,1.8078,1.3037 

B2I  -  2 

1500 

1300 

1260 

0.13 

1 

2.952 

0.4589 

0.4569 

B2I  -  3 

2500 

2122 

2080 

0.15 

1 

23.474 

1 

.3389 

1.3746 

B2I  -  4 

3000 

2504 

2477 

0.16 

1 

— 

— 

0.3732 

2 

47.00 

1 

.8904 

1.8984 

B3II  -  1 

2000 

1733 

1680 

0.13 

1 

6.079 

0 

.6743 

0.6256,0.8153 

B311  -  2 

1500 

1377 

1295 

0.08 

1 

3.745 

0.6992 

0.5377 

B3II  -  3 

2500 

1916 

1987 

0.23 

1 

15.21 

1 

.2054 

0.9113 

2 

9.718 

0.7702 

B3I I  -  4 

3000 

2755 

2590 

0.08 

1 

36.36 

1 

.6971 

0.4215,2.2215,1.3068 

Ball  -  5 

3200 

2929 

2758 

0.08 

1 

15.84 

0.6517 

0.8064 

2 

40.24 

1 

.6555 

1 . 7485 

NOTES:  Min.  Load  = 

0.1  Max. 

Load 

Cn 

(Slope)  (nE4/3Y2/3 )(Bw)2 

APmean  ^initial 

+  APfinal 

)/2 

B 

(AP)2 

=  (0.9Pxi  ♦ 

0.9Pxf)/2 

’  °-45(Pxi 

+  Pxf} 

For  B 3  -  1 : 

Range  Fall 

-Off  Ratio  e  r  -  (*l»1nit1al  - 

c„ 

=  (Slope)(w) ( 28 ) **/ 3 ( 220 ) 2/ 3 ( 0 . 3  X  1.8)2(1010) 

AP 

initial 

B 

( 1 7 1 0 ) 2 

=  Max.  Load  Fall-Off  Ratio 

=  Slope 

'  2839 

X  1010  = 

0.7242 

_  (1710) 

2 

with  data  from 

Table  1. 
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Table  5.  CTS  INTEGRAL  ICTS2n  BASED  ON  ASTM  STAND¬ 
ARD  E399-70T  FORMULA  FOR  R  AND  a  DECREASING  LOAD 
RANGE  aP  =  &PmQa„  (-0.6  a/w  +  1.98).  CONDENSED 
mean  TABLE. 


x  =  a/w  5CTS2D  "  0> 000734 


0.220000 

0 

0.225000 

0.000899 

0.230000 

0.001063 

0.235000 

0.001227 

0.240000 

0.001391 

0.245000 

0.001554 

0.2500p0 

0.001716 

0.255000 

0.001878 

0.260000 

0.002039 

0.265000 

0.002198 

0.270000 

0.002357 

0.275000 

0.002515 

0.280000 

0.002671 

0.285000 

0.002826 

0.290000 

0.002980 

0.295000 

0.003132 

0.300000 

0. 003283 

0.305000 

0.003432 

0.310000 

0.00358D 

0.315000 

'0.003726 

0.320000 

0.003870 

0. 325000 

0.004013 

0.330000 

0.004154 

0.334999 

0.004293 

0. 339999 

0.004431 

0.344999 

0.004566 

0.349999 

0.004700 

0.354999 

0.004832 

0.359999 

0.004963 

0.364999 

0.005091 

0.369999 

0.005218 

0.374999 

0.005343 

0.379999 

0.005466 

0.384999 

0.005588 

0.389999 

0.005708 

0.394999 

0.005826 

0.399999 

0.005943 

0.404999 

0.006058 

0.409999 

0.006171 

0.414999 

0.006283 

0.419999 

0.006393 

0.424999 

0.006502 

0.429999 

0.006608 

0.434999 

0.006714 

0.439999 

0.006818 

0.444999 

0.006920 
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Figure  19.  Theoretical  comparison  of  second-power  crack  growth  curves  corresponding  to  a  constant  range  of  loading 
and  the  gradual  decrease  in  range  of  loading  AP  *  (-0.6  a/w  +  1.198),  of  the  Lal-Weiss  tests. 


FOURTH-POWER  VS  SECOND-POWER  CTS  INTEGRAL  REPRESENTATION 

The  fact  that  the  CTS  integral  vs  N  forms  straight  line  segments  shows  that  the 
growth  rate  does  increase  as  the  square  of  the  crack  tip  stress  intensity  range,  AK. 
As  further  support  of  this,  a  CTS  integral  was  made  up  for  rate  of  growth  per  cycle 
proportional  to  the  fourth  power  of  AK,  i.e.,  m  =  4,  which  has  been  advocated  by 
some  authors  and  is  not  far  from  m  =  3.7  given  the  Lal-Weiss  analysis. 

From  Equations  (2)  and  (25)  it  is 


ICTS4 


/ 

X 

o 


dx 

x2 (29.6  -  185. 5x  +  655. 7x2  -  1017x3  +  638. 9x4)4 


(31) 


Values  of  this  fourth-power  CTS  integral  were  computed  and  tabulated  in  con¬ 
densed  form  in  Table  6.  The  results,  using  this  growth  law,  are  shown  in  Figure  20. 
Here  a  straight  line  was  drawn  to  represent  the  second-power  CTS  integral  vs  cycles 
of  growth.  At  selected  numbers  of  cycles  the  CTS  integral  values  were  read  off  this 
second-power  line.  The  second-power  CTS  integral  tables  were  consulted  to  find 
corresponding  values  of  a/w  and  these  values  were  then  used  to  find  values  of  the 
fourth-power  CTS  integral  from  the  fourth-power  integral  table. 


Table  6.  CTS  INTEGRAL  ICTS4  BASED  ON  ASTM  STANDARD  E399-70T  FORMULA  FOR 
K.  X  =  a/w,  n  =  4.  [SEE  EQUATION  (31).]  CONDENSED  TABLE. 


These  fourth-power  CTS  integral  values  were  then  plotted  on  the  same  graph 
(Figure  20)  as  the  second-power  growth  line  but,  to  make  the  initial  slopes  of  the 
two  curves  the  same  for  comparison,  unit  distance  along  the  fourth-power  CTS  inte¬ 
gral  scale  was  made  26.98  times  as  long  as  on  the  second-power  scale.  This  constant 
is  of  course  the  ratio  of  the  integrands  (i.e.,  dl/dx)  of  I  g£  anc*  ^CTSA  at  t*ie 
lower  limit  of  these  integrals,  x  =  0.2. 

It  is  obvious  from  Figure  20  that  data  plotted  using  a  fourth-power  CTS  inte¬ 
gral  when  the  crack  was  actually  growing  as  a  second  power  would  form  a  curve  with 
its  concave  face  pointing  down,  instead  of  plotting  as  a  straight  line  as  it  would 
if  the  data  followed  a  fourth  power  law. 
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CTS2  a/w  a/w  CTS4 


Figure  20.  Theoretical  comparison  of  second-  and  fourth-power  crack  growth  curves. 

For  the  reverse  case,  using  a  second-power  CTS  integral  when  the  crack  was 
actually  growing  by  a  fourth  power  the  data  would  form  a  curve  with  its  concave  face 
pointing  up.  It  is  only  in  the  case  where  the  CTS  integral  is  of  the  correct  power 
that  the  CTS  integral  vs  N  curve  has  a  straight  line  relationship,  and  on  this 
basis,  the  correct  power,  according  to  the  plots  of  the  experimental  data,  appears 
to  be  two. 


EFFECT  OF  USING  OTHER  FORMULAS  FOR  K  THAN  THE  ASTM  STANDARD  E399-70T 

POWER  SERIES  LAW 

Other  formulae  than  that  in  ASTM  Standard  E399-70T  for  computing  K  for  the  CTS 
have  been  proposed,  such  as  those  in  Table  7.  To  investigate  their  effect  on  the 
matching  of  the  crack  growth  law,  second-power  CTS  integrals  were  made  using  these 
other  relationships  for  AK  (Appendix  A).  A  straight  growth  line  was  drawn  to  repre¬ 
sent  a  hypothetical  CTS  integral  vs  N  line  resulting  from  using  the  ASTM  Standard 
equation  for  AK  in  the  second-power  CTS  integral.  At  selected  numbers  of  cycles  the 
ASTM  Standard  CTS  integral  values  were  read  off  of  the  line.  The  CTS  integral 
tables  were  consulted  to  get  a/w  values,  which  were  then  used  to  find  the  CTS  inte¬ 
gral  values  that  were  computed  as  a  function  of  a/w  using  the  other  AK  equations. 
These  CTS  values  plotted  against  corresponding  values  of  N  show  how  data  giving  a 
straight  line  according  to  the  Standard  formula,  would  look  if  the  substitute 
formulae  for  K  were  used  instead. 

The  resultant  plot  is  shown  in  Figure  21.  The  CTS  integral  scales  for  each 
growth  curve  have  been  displaced  for  clarity.  From  the  plot,  one  can  see  that  there 
is  little  difference,  aside  from  a  slight  curvature,  between  the  shape  of  the  growth 
curves  that  resulted  from  using  different  equations  for  computing  AK  in  the  second- 
power  CTS  integral.  If  the  scatter  in  an  actual  experimental  set  of  data  were  taken 
into  account,  the  differences  might  be  considered  so  negligible  that  the  data  could 


Table  7.  FOUR  SECOND-ORDER  PADE  APPROXIMANT  CTS  FORMULAE  FOR  K  vs  a/w  =  x  REPRESENT¬ 
ING  AND  COMPARED  WITH  ASTM  STANDARD  E399-70T  POWER  SERIES,  TABLES  OF  BOWIE 
AND  SRAWLEY  AND  THE  FORMULAE  OF  BEEUWKES 


Kw1/2  B/P  vs  a/w 


x  =  a/w 

ASTM 

E399-70T 

Pade 

Srawley 
(H/w)  =  0.556 

Pade 

Bowie 

Pade 

Beeuwkes 

Pade 

0.2 

5.19 

5.19 

4.71 

4.71 

4.52 

4.52 

4.35 

4.35 

0.3 

5.85 

5.99 

6.03 

5.92 

5.48 

5.53 

5.57 

5.56 

0.4 

7.32 

7.32 

7.57 

7.57 

7.08 

7.24 

7.24 

0.5 

9.60 

9.59 

9.86 

9.95 

9.56 

9.56 

9.71 

9.73 

0.6 

13.54 

13.54 

13.72 

13.72 

13.70 

13.77 

13.73 

13.73 

0.7 

21.43 

21.58 

21.59 

21.59 

21.24 

21.16 

0.8 

37.42 

37.42 

41.01 

41.03 

— - 

(38.77) 

39.03 

39.03 

Standard 


B  wi/2 

Srawley 

K  »  -- 
B  W1/2 

Bowie 

K  = 

B  w1/2 

Beeuwkes 

K  *-P 
B  w1/2 


4,5170  -  4.827Sx  +  7.0528x2 
1  -  1.3867X  +  0.3867x2 


3.0532  -  0.0167X  -  3.7528x2 
1  -  2.1528X  +  1 . 1 528x2 


3.4152  -  1.6577X  +  3.6408x2 
1  -  1.5736X  +  0.5376X2 


2.7168  +  0.8872x  -  3.0716x2  ' 
1  -  2. 01 61  x  +  1.0161x2 


still  be  taken  to  correspond  to  a  straight  line  whichever  formulae  might  be  used. 
In  other  words,  the  different  formula  proposed  for  computing  AK  may  not  correspond 
to  materially  different  growth  curves  in  the  range  of  a/w  considered,  insofar  as 
comparison  with  experimental  data  scatter  or  variability  is  concerned.  They  do 
suggest  however,  that  in  our  analysis  of  a  crack  growth  experiment,  the  apparent 
occurrence  of  a  locus  of  crack  growth  made  up  of  lines  having  very  small  slope  dif 
ference  should  be  disregarded  in  favor  of  a  single  straight  line  since  a  general 
law  of  growth  must  be  inferred  from,  or  compared  with,  nonideal  experimental  data. 
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Figure  21 .  Theoretical  comparison  of  second-power  crack  growth  curves  derived  from  stress  intensity  K  formulae  for 

Beeuwkes,  Bowie,  Srawley,  and  ASTM  Standard  E399-70T. 


APPARENT  DISCREPANCY  m  =  2  vs  m  =  3.7  RESOLVED;  da/dN  vs  AK  PLOTS 

Our  analysis  has  demonstrated  that  the  crack  growth  rate  of  Trip  Steels  can  be 
completely  represented  by  second-power  relations  despite  the  fact  that  Lal-Weiss, 
using  accepted  experimental  procedures,  obtained  a  markedly  higher  power,  m  =  3.7. 

The  explanation  of  this  may  be  found  in  log  da/dN  vs  log  AK  plots.  Figures  22, 
23,  and  24,  of  all  the  data  on  a  material,  in  which  the  Lal-Weiss  points,  that 
correspond  to  grown  crack  lengths  of  only  0.2  in.,  are  encircled.  In  these  plots 
it  is  seen  that  the  growth  rates  in  all  experiments  were  proportional  to  (AK)2,  even 
though  higher  growth  rate  curves  corresponded  to  larger  loads.  The  Lal-Weiss  data 
points  as  a  group  do  not  fit  any  experimental  growth  curve;  the  points  jump  from 
one  curve  to  another  but,  since  all  correspond  to  a  fixed  crack  length,  their 
locus  does  have  the  same  slope  as  the  envelope  of  all  data  points.  Except  for 
Specimen  Bj  #1,  Heat  2321,  Type  I,  practically  all  data  points  lie  within  an 
envelope  consisting  of  two  parallel  lines  having  a  slope  of  four  and  separated  from 
each  other  in  accordance  with  PRESENT  TYPE  OF  ANALYSIS. 

The  AK  and  da/dN  values,  Table  8,  on  these  figures  were  computed  from  crack 
growths  separated  by  0.05-in.  increments  taken  from  smooth  curves  drawn  through  the 
experimental  data  as  discussed  following  Equation  (28)  and  shown  on  the  sample  data 
sheets  from  Lal-Weiss  (Figure  2).  They  were  computed  for  m  ■  2,  as  follows: 
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Table  8.  da/dN  AND  aK  FOR  FIGURES  22,  23,  AND  24 


Slope 

at  *  '°6 

AK 

Slope 

Ix  106 

aK 

Thk 

a, in. 

Specimen 

xl  0” 

in./cy. 

ksi  -/in. 

Specimen 

xl0° 

in./cy. 

ksi  ./in. 

0.40 

4.88 

2.43 

22.3 

3.66 

1.82 

16.5 

0.425 

2.49 

22.6 

1.87 

16.7 

0.45 

2.57 

23 

1.93 

16.9 

C 

0.50 

O 

9.24 

5.26 

23.9 

O 

2.08 

17.6 

cm  m 

o 

o 

0.55 

Cv. 

5.81 

25.1 

CM 

2.30 

18.5 

o 

0.60 

•— <  II 

6.52 

26.6 

►— <  II 

2.58 

19.6 

CO  CL 

CO  CL 

II 

0.65 

<3 

6.07 

4.85 

28.3 

2.29 

2.93 

20.9 

CO 

0.70 

5.55 

30.3 

2.09 

22.3 

0.75 

6.38 

32.5 

2.41 

23.9 

0.80 

7.37 

34.9 

2.78 

25.7 

0.40 

6.59 

3.28 

23.2 

i 

2.94 

1.46 

17.3 

0.425 

3.36 

23.5 

1.50 

17.5 

0.45 

20.8 

10.9 

23.8 

1.55 

17.8 

c 

0.50 

m 
«—  CO 

11.9 

24.8 

O 

1.67 

18.5 

LO 

CO 

i O 

CM  <0 

0.55 

' 

13.1 

26.1 

CM 

•  •  >— 

1.85 

19.5 

d 

0.60 

CM 

CO  O- 

14.7 

27.6 

•— «  II 

CN 

2.07 

20.6 

< 

CO  O- 

0.65 

15 

12.0 

29.4 

<J 

2.35 

22 

0.70 

13.7 

31.4 

2.69 

23.5 

0.75 

15.8 

33.7 

3.09 

25.2 

0.80 

18.2 

36.2 

3.57 

27.1 

0.40 

7.35 

3.66 

20.9 

2.88 

1.43 

15.7 

0.425 

3.75 

21.2 

1.47 

16.3 

0.45 

3.87 

21.5 

1.51 

16.6 

c 

0.50 

r-v. 

4.19 

22.4 

1.64 

17.2 

m 

CO 

0.55 

d 

i  r^. 
to 

4.62 

23.5 

l  CT> 

C\J 

1.81 

18.1 

O 

0.60 

5.64 

5.19 

24.9 

2.03 

19.2 

ii 

CO 

0.65 

CO 

O- 

< 

4.51 

26.5 

CO 

CL 

<3 

2.30 

20.4 

0.70 

5.15 

28.3 

2.63 

21.8 

0.75 

5.92 

30.4 

3.03 

23.4 

0.80 

6.85 

32.7 

3.50 

25.2 

For  B  3 1 1-2  and  B  3 1 1 - 3  Lal-Weiss  GCL  vs  N  curve  was  slightly  modified  to  make 
*CTS2  vs  N  a  slight  line. 

AK  =  (aP/B  w‘/2)(x1/2)(29.6  -  185. 5x  +  655. 7x2  -  1017x3  +  638. 9x** )  from  ASTM  E399-70T 
x  5  a/w  (Figure  1 ) 

S  =  W  (--p7^wi/2)2  X  S1°Pe  Whpre  S,0pe  =  (6lCTS2/aN>  on  Lal-Weiss  Smooth  Curves  of 
Grown  Crack  Length  ( GLC )  vs  Cycles  (N)  data 


0.315  in.  B  =  0.2B5  in.  B  =  0.100  tn 


Table  8  (Cont).  da/dN  AND  AK  FOR  FIGURES  22,  23,  AND  24 


Slope 

—  x  106 
dN  x  ,u 

aK 

Slope 

dN  x  10& 

AK 

a,  in. 

Specimen 

xlO8 

in./cy. 

ksi  /in. 

Specimen 

xlO8 

in./cy. 

ksi  ./in. 

0.40 

6.07 

28.6 

10.3 

5.13 

34 

0.425 

12.2 

6.22 

29 

5.25 

34.4 

0.45 

6.42 

29.4 

5.42 

35 

0.50 

CM 

6.95 

30.6 

CM 

25.6 

14.6 

36.4 

cn  on 

O 

0.55 

1  CM 

18.9 

11.9 

32.2 

VO 

1  CM 

16.1 

38.2 

0.60 

>— •  II 

13.3 

34.1 

►—*  II 

18.1 

40.5 

cD  a. 

CO  Q_ 

0.65 

25.3 

15.1 

36.3 

< 

20.5 

43.1 

0.70 

23.1 

38.8 

23.4 

46.1 
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For  B3II-2  and  B  3 1 1 - 3  Lal-Weiss  GCL  vs  N  curve  was  slightly  modified  to  make 
ICTS2  vs  N  8  straight  line. 

AK  =  (AP/B  w‘/2)(x1/2){29.6  -  185. 5x  +  655. 7x2  -  1 01 7x3  +  638. 9x4)  from  ASTM  E399-70T 
x  =  a/w  (Figure  1 ) 

^  S  W  (  aP/B^W2)2  *  S1°Pe  Wh8re  S1°Pe  =  !fllCTS2/AN)  on  Lal-Weiss  Smooth  Curves  of 
Grown  Crack  Length  (GLC)  vs  Cycles  (N)  data 
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Table  8  (Cont).  da/dN  AND  AK  FOR  FIGURES  22,  23,  AND  24 
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l.e. , 


ICTS2  |a  "  ICTS2  la 
_ o 

N  -  N 

o 


(6) 


where  la  means  at  a. 

The  da/dN  values  so  computed  for  a  =  0.4  +  0.2  in.  (i.e.,  grown  crack  length  = 
0.2  in.)  were  in  substantial  agreement  with  the  values  obtained  by  Lal-Weiss  so  that 
no  graphical  distinction  between  the  sources  of  da/dN  data  is  called  for. 

However,  our  AK  values  which  ar^  computed  from  ASTM  Standard  E399-70T,  as  are 
the  Lal-Weiss  values,  are  about  1.25  times  as  great  as  the  Lal-Weiss  values.  Their 
f  (a/w)  =  5  in  AK  = (AP/B^/w)f (a/w)  of  this  standard  according  to  Figure  2(a), 
whereas  f(a/w)  =  6.26  actually,  for  a/w  =  (0.4  +  0.2)/1.8  where  0.4  is  the  initial 
crack  length  to  which  the  "grown"  crack  length  of  0.2  in.  is  to  be  added  to  get  the 
total  distance  from  the  action  line  of  the  pins  of  the  loading  fixture  to  the  crack 
front.  Thus,  our  values  of  AK  were  used  in  the  plots  (Figures  22,  23,  and  24). 

(Use  of  the  Lal-Weiss  AK  values  would  simply  shift  all  points  to  the  left  a  constant 
amount  -  an  inconsequential  difference  as  far  as  the  Lal-Weiss  determination  of  m 
is  concerned. ) 


CONCLUSIONS 

Because  of  discontinuities  in  da/dN  that  are  not  readily  observable  on  crack 
growth  curves  drawn,  as  is  customary,  to  Idealize  periodic  observations,  it  is  best 
to  match  the  actual  undifferentiated  data  with  a  function  (possibly  segmented)  hav¬ 
ing  plausible  characteristics  and  determine  da/dN  from  this  function. 

The  crack  growth  curves  of  the  Trip  Steels  analyzed  here  closely  follow  the 
Beeuwkes  formulation.  PRESENT  TYPE  OF  ANALYSIS,  which  has  been  used  successfully 
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on  many  materials  and  according  to  which  there  are  discontinuities  in  da/dN  and 
except,  at  these,  da/dN  is  proportional  to  (AK)2.  This  agreement  follows  from  the 
experimental  results  whether  or  not  there  may  have  been  nonuniform  material  or  a 
crack  front  that  changes  shape.  Appendix  A  is  supplied  herewith  to  facilitate 
plotting  growth  curves  by  this  formulation. 

The  apparent  discrepancy  between  the  Lal-Weiss  relation  da/dN  ~  (AK)3'7  and  our 
da/dN  ~  Const,  x  (AK)2  is  that  their  relation  does  not  correspond  to  a  rate  of 
growth  along  an  actual  growth  curve  while  ours  does.  da/dN,  AK  points  that  Lai  and 
Weiss  plotted  to  get  their  relation,  jump  from  a  point  on  one  growth  curve  to  a 
corresponding  point  on  another,  and  so  on,  and  the  change  in  rate  for  a  given  range 
of  AK,  along  such  a  relation  is  not  the  same  as  the  change  on  a  single  growth  curve. 
Their  m  =~3.7  (~4)  is  the  slope  of  a  line  connecting  corresponding  points  on  the 
different  segments  (all  having  a  slope  of  m  =  2)  of  the  log  da/dN  vs  log  AK  diagram 
and  is  parallel  to  the  envelopes  connecting  line  segment  ends  in  this  diagram  since 
all  the  segments  seem  to  have  the  same  length. 

Although  there  is  no  growth  curve  for  which  da/dN  is  proportional  to  (AK)1*,  one 
may  argue  (in  conformance  with  the  Beeuwkes  formulation)  that  in  a  certain  statis¬ 
tical  sense  the  fourth  power  law  may  be  considered  correct.  Namely,  for  the  average 
performance  of  a  large  number  of  samples  of  a  material  or  in  the  case  that  the 
growth  or  increase  in  AK  to  double  AK  renders  AK  still  so  small  as  to  be  inconse¬ 
quential  for  some  application;  under  constant  stress  this  means  that  the  growth  in 
a  crack  length  to  double  that  length  makes  the  length  still  too  small  to  be  con¬ 
sidered  important.  However,  even  if  the  length  is  not  negligible,  there  is  a  fourth 
power  envelope  line  to  the  growth  curves  that  may  confidently  be  arbitrarily 
assumed  for  conservative  design  to  be  the  growth  curve.  (See  PRESENT  TYPE  OF 
ANALYSIS  and  SECOND-POWER  CTS  INTEGRAL  AND  ANALYSIS;  AK  BY  ASTM  STANDARD  E399-70T.) 
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APPENDIX  A.  SECOND-POWER  CTS  INTEGRALS 


CTS  Integrals,  IqtS2  vs  x  -  a/w*  f°r  Constant  Loading  Range  for  ASTM  Standard 
E399-70T  Power  Series  for  K,  Second-Degree  Pade  Formulae  for  Beeuwkes  and  Bowie  K 
and  Third-Degree  Pade  Formula  for  Srawley  K,  follow  in  Table  A-l. 
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